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Abstract 

Specific dark energy models, leading to the Little Rip (LR) cosmol- 
ogy in the far future, are investigated. Conditions for the occurrence 
^ ■ of LR in terms of the parameters present in the proposed equation of 

bi) state for the dark energy cosmic fluid are studied. Estimates about 

the time needed before the occurrence of the rip singularity in the 
standard LR model and the model in which the universe approaches 
the de Sitter space-time asymptotically, are given. 



1 Introduction 

The discovery of the accelerating universe has led to the appearance of new 
theoretical models (for recent reviews, see [HE])- The cosmic acceleration 
can be explained via the introduction of dark energy (recent references are 
[3111]), with strange properties like negative pressure and/or negative entropy. 
According to the latest supernovae observations the dark energy accounts for 
73% of the total mass energy of the universe (see, for instance, [5]). 

One characteristic feature of the new cosmology is the Big Rip (BR) 
phenomenon, i.e. a singularity of the universe to be encountered in a finite 
time. Pioneering papers on the Big Rip were given in P [3 El El ED]- The 
Big Rip means that one or more of the physical quantities go to infinity in 
a finite time t in the future. Mathematically, it implies divergent integrals 
following from the Friedmann equations. 
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A softer variant of the future singularity phenomenon is the so-called 
Little Rip (LR), characterized by an energy density p increasing with time 
but in an asymptotic sense, so that an infinite time is required to reach 
the singularity. The LR scenario was proposed in It corresponds to 

an equation of state parameter w < —1, but — )■ — 1 asymptotically. A 
viscous generalization of the LR theory was recently given in [12]. The energy 
density grows with time but not fast enough for the occurrence of the Big 
Rip singularity. The LR models describe intermediate evolutions between an 
asymptotic de Sitter expansion and a BR evolution. 

The so-called Pseudo-Rip (PR), proposed in [13], is a third variant of 
the theory. This interesting possible scenario is related to the LR cosmology 
when the Hubble parameter tends to infinity in the remote future [T^ 



A model discussing LR in multivariate Gauss-Bonnet theory with dilaton 
was constructed in [19]. The LR cosmology may be described in terms of a 
general fluid with a complicated equation of state [3 [20] . The nonsingular 
LR cosmology may lead to structure disintegration in the future (galaxies, 
the Sun system, etc), similarly as the BR. 

In the present paper we study the influence of the time dependent pa- 
rameters w and A in the equation of state (see equation ([5]) below) upon the 
occurrence of LR/PR in various cosmological models. Estimates for the time 
required for disintegration of gravitational bound systems are given. 

2 Equation of state with time-dependent w 
and A in the Little Rip model 

Let us consider the following LR model with a given Hubble parameter H: 



Here it is natural to associate t = with the present time, so that Hq becomes 
the present-time Hubble parameter. The Friedmann equation for a spatially 
flat universe is 



[161 in [18], 



H{t) Hoo < oo, t oo. 




Hit) = Hoe^\ Ho>0, A > 0. 



(2) 




(3) 
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where p is the energy density, H = a/a the Hubble parameter, a{t) the scale 
factor, and = SnG with G being Newton's constant. Thus 



P=^ff^ (4) 



We assume now that our universe is filled with an ideal fiuid (dark energy) 
obeying an inhomogeneous equation of state [21] 



p = wit)p + A{t), (5) 

where p is the pressure and w{t),A{t) are time-dependent parameters. This 
form is related to the one used in our earlier paper [22] ; cf . also [23] . 
The energy conservation law is 

p + 3H{p + p) = 0. (6) 

Taking into account (j4l)-(l6]) we obtain 

3 



^H' + 3H 



\l + w{t)]H' + A{t) 



(7) 



Solving equation ([7]) with respect to w{t) we have 



2 A A(t)K2 

Let us investigate various cases: 

1. Assume that the parameter A does not dependent on time, A(t) = Aq. 
Then w{t) —1 asymptotically from below, in the far future. Thus we see 
that for an ideal fiuid obeying equations ([5]) and ([8]) the LR scenario is found. 

If instead the parameter w does not depend on time, w{t) = wo, we can 
solve equation (JTj) with respect to Ait) to get 

A{t) = -^[2X + 3Hil + wo)]. (9) 

As wq < —1 in order to correspond to a dark fiuid, we see that A(t) — )• oo in 
the far future. The LR behavior in this case is caused by the parameter wo- 
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2. Let us now consider instead of ([2]) another LR model [T2] where the 
Hubble parameter is given as 

H{t)=Hoe'^'"\ (10) 

Here Hq, C and A are positive constants. (Note that with this definition, 
the present-time value of the Hubble constant becomes i^(0) = Hqc^.) In 
the paper [TT] it was shown that this model may be a realistic one, as it is 
compatible with the observations. 
We now find 

P - -^^0 e , 

and in view of ( |T0|) and ( ITTi) equation ([6]) becomes 

3 



Ho 



K 



^[l+wit)]H' + A{t) 



0. (12) 



Let us suppose that the parameter A{t) is proportional to the square of the 
Hubble parameter [21], that is 

m = iH\ (13) 

with 7 a constant. Taking ( IT3!) into account, we solve equation (fT2|) with 
respect to w{t) to obtain 

Now writing the parameter w{t) in the form 
with 5 a positive constant, we obtain from f fT2|) 

As shown in ^llj we can easily find models in which there are more compli- 
cated behaviors of H, such as 

H = Hq exp (Co exp Ci (exp C2 (exp ... (C„ exp(At))))) , (17) 
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where Cq, Ci, are positive constants. In this case the energy conserva- 
tion law takes the form 



Ho Ho 

where we have defined Ln as 



■[l+w{t)]H^ + A{t) 



0, 



Ln- 



H 
^Ifo 



In 



— In — 

Co Ho 



In 



-...hi 



k-l 



— In — 

Co Ho 



with k G (1, n) and n & N . 

Using ( !T3|l to solve equation (fT8|) with respect to we obtain 



2\ H ^ H 

w{t) = In — Ln K 

3 H Ho Ho 3 
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1. 



Analogously, by generalizing f ITB]) we obtain 

2X^^H^H 35 
At = -— ^In— Ln— + -. 



(18) 



(19) 



(20) 



3. Now go on to consider a Pseudo-Rip model [13J with a different behavior 
oiH [}A\: 

H{t) = Ho- Hie-^', (21) 

where Ho, Hi and A are positive constants. We assume that Ho > Hi when 
t > 0. Since the second term decreases with increasing t, the universe ap- 
proaches asymptotically the de Sitter space-time with Hubble constant Ho 
[Ho here means the same as H^o in ([1])). 
We now find 



P 



iHo-Hie 



-Xt\2 



P 



6A 



H{H - Ho 



Using f|T8|) we can rewrite the energy conservation equation as 



6A 



H{H - Ho) + 3H 



:[l+w{t)]H^+A{t) 



0. 



(22) 



(23) 



We will now investigate this kind of PR, in analogy with the earlier models 
P and (HDD. 
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First, if we assume A{t) = Aq, then we find 

wfiich shows that the PR behavior is determined by the parameter Aq. 
Next, if we take A(t) = 'jH'^, we find 

-W = -3A^^- — -1, (25) 

showing that the PR is connected with the Hubble parameter (!2T|) . 

Solving A(t) from (125]) and taking w(t) = Wq, we see that the PR is 
determined by the parameter wq: 

A(t) = - — {H - Ho) - 4(1 + ^o)^'. (26) 
In view of (|T5|) this means 

9 A S 
m = --{H - H,) - -5. (27) 

Thus, we have presented the appearance of LR and PR from the equation of 
state dSl). 



3 The inertial force in Little Rip cosmology 

During the universe expansion the relative acceleration between two points 
whose distance is / is equal to /a/a, where a is the scale factor. A particle 
with mass m at a given point will be subject to an inertial force [13] 

Fi^ = ml- = ml{H + H"^). (28) 
ct 

Let us assume that two particles are bound by a constant gravitational force 
Fq. If Fin > and Fi^ > Fq, the particles become unbound. This is the 
Big Rip phenomenon caused by the accelerating expansion; the Sun system 
or the galaxies in the universe may become gravitationally unbound. It is 
convenient to define the dimensionless parameter [131 [TH] 

2p(a) + p'(a)a a 

Fin = ^ ^ ^ ^ = 6 , (29) 

Po apo 
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where po is the dark energy density at present. 

Now consider the LR model as determined by ([2]). The inertial force is 



Fi^ = ml{XHoe^' + Hy^'). (30) 

Hence, at time t — )■ +00 the inertial force Fin — ^ +00. This characterizes 
the LR; under certain conditions a disintegration may occur for gravitational 
force structures. 

Let us investigate the influence from the parameters w and A in the 
equation of state for dark energy upon the time tLR needed for disintegration 
in the LR model governed by ([2]). In the equation of state ([5]) we take 

w(t) = Wo, 

Wo < -1, (31) 
A{t) = Ao, 

where Wo and Aq are constants. 
Assume now that 

A' 

The energy conservation law becomes 



^0 = -1 - TT^. (32) 



(XHofe^^' - 2Ao/s:'Aifoe^* - = 0- (33) 

From this we find the disintegration time tiR'- 

1 (1 + V2)AoK^ 
^^^=A^" Wo ' 

Parameter estimates show that the Rip takes place after some billions of 
years, and that the Sun-Earth system disintegrates when Fin is of order 10^^ 
years [15] . 

As another example, consider the PR model (1T71) . The inertial force (HM 
is equal to [13] 

Fin = ml [Aifie-^* + {Hq - H^e-^'f] . (35) 
In this case the inertial force is limited: Fin ~^ mlHl when t — )■ +00. 
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As the time t increases, the universe develops to an expanding de Sitter 
space-time determined by the parameter Aq. Choose the relationship between 
Wo and Aq in equation ([5]) in the form 

^o = -l-^. (36) 
Then the energy conservation law is written as 

.^Ao^e-- + 2(A-^)=0, (37) 
and one finds for the Rip time tpR 

It is thus possible, from the relations (IMj) and (EH]), to estimate the span of 
time tiR or tpn needed before the system becomes gravitationally unbound. 



4 Conclusions 

In summary, we have presented dark energy models with an inhomogeneous 
equation of state, cf. equation ([5]) above, in which Little Rip or Pseudo- 
Rip behavior is encountered in the far future. We have shown that the LR 
cosmology is caused exponentially, determined by the parameters A or w. 
It is of interest to note that the disintegration of bound structures in the 
LR/PR models may occur for physically acceptable choices of parameters. 
We have given estimates for the disintegration times, in case of the standard 
model for LR, and also for the case of the asymptotic de Sitter universe. 

As a final comment, we mention that it is possible to generalize the LR 
theory so as to take into account the bulk viscosity in the dark fiuid. Such 
a formulation was developed in [?T2j. One convenient choice for the equation 
of state, which permits an analytic solution to be found, is the following, 

p = -p-A^-^o. (39) 

Here A is a positive constant (with dimension cm~^ in geometric units), 
and the last term C,o is a constant representing the infiuence from the bulk 
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viscosity. Some calculation (details omitted here) shows that p(t) increases 
exponentially for large values of t. This is the characteristic feature for the 
LR again, now with the inclusion of viscosity. 



Acknowledgement. We thank Professor Sergei Odintsov for very useful 
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